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We study the recent PHOBOS data on the pseudorapidity density of inclusive 
charged particles in centrality-binned d+Au collisions at y/s^N = 200 GeV. It 
£ — ■ appears that one can understand the increasing forward-backward asymmetry in 

the data by assuming that the entire distribution "shifts" backwards in rapidity 
according to the initial-state kinematics, while the total multiplicity scales linearly 
with Np ar t/2. Two models are explored, both of which achieve a reasonable descrip- 
tion of the available data for rapidities sufficiently far from the projectiles (\n\ < 3 
at the top RHIC energy). One model uses PYTHIA as the underlying distribution, 
to allow a straightforward mapping of a rapidity shift back to pseudorapidity space. 
The other model is a simple analytic calculation based on functions inspired by 
Landau's hydrodynamical model. The apparent success of these simple pictures to 
describe the bulk of particle production over \n\ < 3 suggests that collective effects 
may be present even in the "small" systems created in p+p and d+Au reactions, ac- 
tive over the full rapidity range. The relationship between this and other theoretical 
approaches are discussed. 



Recent PHOBOS data [1] have shown the pseudorapidity distribution of charged 
particles for minimum-bias d+Au collisions at the full RHIC energy snn = 
200 GeV as well as the centrality dependence of the distributions. In this work, 
we will show that the bulk particle production in p+A and d+Au collisions 
over a range of energies can be described by a very simple model with only 
two assumptions: 

• The total multiplicity in p+A and d+Au collisions is linear with the num- 
ber of participant pairs (N part /2) relative to p+p collisions (similar to the 
"wounded-nucleon model" [2,3]. 

• The angular distribution in p+A and d+Au collisions are approximately the 
same as p+p collisions (if one avoids the very far-forward rapidity regions) 
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Fig. 1. The center-of-mass rapidity shift Ay vs. v, the ratio of number of participants 
in the nucleus to that in the smaller projectile. The same function shifted up by 0.3 
units is also shown,to account for "extra" rapidity shifts due to the geometry or to 
transverse dynamics. 

but in a "shifted" reference frame given by the center-of-mass system denned 
by the participants in the smaller and larger projectiles. 

Simple kinematics gives the result that the rapidity of the center-of-mass (CM) 
frame in a p+A collision, where the proton interacts with a "tube" of v nu- 
cleons in the nucleus, is 



In a d+Au collision, one can approximate the geometry to be the independent 
collision of the proton and neutron with half of the participants in the gold 
nucleus, or a full collision of the 2 nucleons in the deuteron with all of the 
participants in the nucleus. In either case, 



In general, these two definitions lead to equivalent physics conclusions, with 



According to the above assumptions, the produced particles should emerge 
with a rapidity distribution similar to the underlying p+p distribution, but 
shifted by Ay relative to the original p+p CM frame. However, it should 
be noted that if the relevant physics involves rescaling xp = 2p T /y/s, where 
dN/dxF is universal with energy for a given initial state geometry, then there 
might be an additional modification to the rapidity distribution, due to the 
transverse momentum of the emitted particles. This is easily seen by the ap- 
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Fig. 2. PHOBOS dN ch /dri distribution for d+Au collisions compared with PYTHIA 
pp simulations scaled up by N part /2 and with the rapidity of each particle shifted 
by Ay =0.7 and 1. 

proximate expression relating y and Xf (good for y away from 0): 



y = y b + ln(x F ) - 




Without a full knowledge of the longitudinal and transverse distributions, it 
is difficult to estimate the complete effect of the dynamical evolution on the 
final state rapidity distributions. One might also consider some additional 
effects from reinter act ions with the spectator matter from the nucleus, sup- 
plying the system with some extra longitudinal momentum in the backwards 
(A) hemisphere. Despite these uncertainties, it will be assumed that any ad- 
ditional modifications to the longitudindal shift (e.g. from effects associated 
with transverse dynamics or spectator interactions), can be incorporated into 
a constant extra contribution to Ay, hereby labeled Sy. The basic relationship 
between the initial-state geometry and the rapidity shift is shown in Fig.l, 
with and without an extra shift of 5y = 0.3. 

The simplest application of the approach described here can be performed 
with the PHOBOS minimum-bias d+Au data, with v estimated to be approx- 
imately 4. This data, shown in Fig. 2, displays a clear rapidity asymmetry, 
with a large peak seen near rj = —2 and an intriguing "plateau" near r] — 1. 
These features are often interpreted in terms of the independent fragmenta- 
tion of the deuteron and nucleus in the forward and backward rapidity regions 
respectively. However, we can use the procedure outlined above to construct 
d+Au distributions from p+p, using PYTHIA 6.161 [4]. First one scales up 
the PYTHIA p+p dN c h/drj distribution up by a factor of N part /2 ~ 4, indi- 



cated by a dotted line. This gives a distribution with similar integral to the 
PHOBOS data, but one which is clearly higher in the forward region and lower 
in the backward region. Then we apply a rapidity shift Ay to each particle in 
the PYTHIA simulation and then recalculate dN c h/dr] fory^y — Ay. Using 
the straightforward formula, Eq. , we would get Ay = 0.7, shown by the dot- 
dashed line. Already we see that this result describes the basic features of the 
PHOBOS data, including the peak at rj = —2, which is now seen clearly as 
a kinematic effect resulting from using pseudorapidity (if) instead of rapidity 
(y). However, it is also shown that assuming 5y = 0.3, giving a total Ay = 1.0, 
gives a somewhat better fit with the data over the measured range. 

This exercise shows that the simple model proposed above can in fact achieve 
a satisfactory quantitative agreement with the data. However, a slight modifi- 
cation (8y = 0.3) can improve this agreement substantially. The next question 
then how far one can test these ideas by varying the collision centrality, which 
changes v. Fortunately, the PHOBOS collaboration has also published data 
comparing d+Au data binned in centrality bins (5 bins, each containing 20% 
of the total inelastic d+Au cross section) [5] , and presented as a ratio compared 
to p+p data from UA5, i.e. the variable 

_ dN d+Au /drj 
v dN p+p /dr) ' 

The BRAHMS collaboration has presented d+Au data binned in centrality 
bins [6], but their centrality method (using particles emitted into \r]\ < 2.2) 
creates severe biases on the final distributions. This makes direct comparisons 
with the models shown here quite difficult without substantial simulation work 
to remove them. 

We have made theoretical calculations of R v to compare with the PHOBOS 
data in two different ways. One is to use the prescription above (N part /2 scaling 
and a Ay shift) on PYTHIA calculations, which has the advantage of allowing 
a full calculation involving the translation from y to r], to see the effect of the 
Jacobian. In Fig. 3, we show these PYTHIA calculations as solid histograms, 
which describe the data reasonably well within \r)\ < 3. The two panels both 
illustrate that the Jacobian essentially cancels in the ratio of d+Au and p+p. 
However, in the comparison of 5y = (which is a parameter-free calculation) 
and 5y = 0.3, it is observed that the non-zero 5y gives a somewhat better 
description of the data. These results thus support and extend the one shown 
in Fig. 2. 

The other model is an analytic calculation based on Landau's hydrodynamics, 
as presented by Carruthers and others [7]. In this approach, the rapidity dis- 
tribution in p+p is a Gaussian with width a 2 = (1/2) ln(s/4m 2 ), a function 
particularly simple to scale by N part /2 and shift by Ay. One can then calculate 
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Fig. 3. The ratio of d+Au and p+p pseudorapidity distributions is shown for a 
range of centralities in d+Au (corresponding to Ref. [5]). The data is compared to 
two calculations, one using PYTHIA (dotted histograms) and the using a simple 
analytic model inspired by Landau's hydrodynamical model (solid line). The left 
and right panels are for Sy = 0.0 and Sy = 0.3, respectively. 

R y , which is the same as but in rapidity space: 



_ Jy parte -(yAy)/* 2 

2 

which gives a exponential dependence of R y with y. For simplicity, we will 
approximate y in the final expression by 77, which appears justified by the 
PYTHIA results, as the r]-y transformation does not seem to have a great 
effect on R n . These results are shown in Fig. 3 as solid lines, and also give an 
adequate representation of the data, within \q\ < 3. 

To test the applicability of the basic approach outlined here to lower-energy 
reactions, the procedure outlined above has been applied to the NA5 data 
set of dN/dy for inclusive produced charged particles measured in 200 GeV 
protons incident on hydrogen, argon and xenon targets [8]. Varying the target 
nucleus varies the number of collisions and allows us to test the predicted 
formula for the rapidity shift, as well as N part scaling. The first panel of Fig. 4 
shows the raw data on dN/dy, the rapidity distribution of inclusive charged 
particles assuming the pion mass. The second panel shows these distributions 
divided by N part /2, where N part — v + 1 and the values of v are taken from 
Ref. [8]. The third panel shows the scaled distributions shifted to the right 
using the formula derived above (with Sy = 0.3). After these transformations, 
a reasonable agreement of the distributions is observed over a large rapidity 
range. It is also observed that the p+A collisions are quite symmetrical over 
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Fig. 4. left) NA5 data on dN/dy for inclusive produced charged particles, shown for 
three different targets, center) The same data scaled by the number of participant 
pairs, right) The same data as a function of y + Ay, with Sy = 0.3 (except for p+p). 

the full range relative to the centroid, even when shifted by up to a unit in 
rapidity. 

Rapidity shifts in the context of "fireball" desriptions of heavy ion collisions 
have been discussed since the late 1970's[9]. However, as far as the author 
is aware, the shifted rapidity distribution for p+A collisions has only been 
incorporated explicitly into one concrete physical model, the Collective Tube 
Model (CTM) [10]. The CTM modelled p+A collisions as being the collision 
of a proton with another particle of the same size as a proton, but with a 
mass and momentum v times larger. This also predicted a shift of all produced 
particles by ln^ 1 / 6 ) (compatible with the simple calculation shown here given 
that v oc A 1 / 3 ). However, to rigorously conserve longitudinal momentum as 
well as total energy, this model also predicted substantially fewer particles 
than observed in a large number of experiments, and a centrality dependence 
(oc Npar t ) which disagrees with the A^^-scaling observed in all existing p+A 
data. 



Shifted rapidity distributions have also been discussed in the context of Lan- 
dau's hydrodynamical model [11]. In Landau's approach, the incoming pro- 
jectiles thermalize in a volume reduced by the Lorentz contraction, which in- 
creases with energy as fast as the CM energy itself. The angular distributions 
are then the consequence of a rapid longitudinal expansion, which can be cal- 
culated in the context of relativistic hydrodynamics. Landau and Belenkij [12] 
performed a detailed calculation on particle production in p+A collisions, 
and found that under reasonable conditions (is < 4), that the total multiplic- 
ity should scale with N part and that the angular distributions should differ 
"only slightly" from that found in nucleon-nucleon collisions. Carruthers [7] 
suggested that asymmetric collisions should lead to "shifted gaussians", but 
detailed phenomenological studies were not made. The salient point to distin- 
guish the Landau approach from the CTM approach is that in the former, one 
would not assume all of the energy in the tube of v nucleons to be compressed 
into a volume the size of a (Lorentz-contracted) nucleon. Rather, it has a lon- 



gitudinal extent approximately v times a nucleoli, such that it's thermalized 
energy density is not dissimilar from a p+p collision. 

The most important issue is that one should not expect the full inclusive 
distribution to be shifted by Ay, as this would manifestly violate energy- 
momentum conservation if extended all the way to beam rapidities, not in 
the least when one shifts the leading particles on the nucleus side to a larger 
rapidity than they had before the interaction. Thus, it is perhaps not surprising 
that the two models presented here fail at large rapidities on both the projectile 
and target sides, as seen in Fig. 3. Of course, the details of these deviations 
may well give useful information either about the thermalization process, or 
even about the particular issues involved with hydrodynamic evolution near 
to kinematic boundaries. 

Finally, we compare the model sketched here with several existing in the lit- 
erature. In effect, we present an alternative explanation to what had been 
suggested by Brodsky et al. [13] to be the result of the multiple-collision pro- 
cess and overlapping "strings" between the colliding nucleons. They predicted 
a linear decrease of R v with r\, while the models here predict an approxi- 
mately exponential decrease. There is also a substantive difference between 
this approach and that advocated recently by Bialas and Czyz, motivated 
by a successful fit of the same data with an updated "wounded nucleon" 
model [14]. In particular, the physical picture advocated here suggests that 
it is the collective system which expands in a shifted frame to make the final 
state rapidity distribution, rather than the independent fragmentation of the 
incoming nucleons. This avoids the situation, detailed by Bialas and Czyz, that 
the forward-going nucleons can fragment into fast-moving particles heading in 
the opposite direction. 

Of course, the limited precision of both the data and the theory precludes any 
decisive choice of any of these models at this moment. The physical picture 
outlined in this work has some appeal having essentially no free parameters 
except the residual shift (Sy), and requires no additional factorization of p+p 
collisions into two fragmentation regions, as suggested by "wounded nucleon" 
approaches. However, it should be noted that explaining this "shift" in the 
overall rapidity distribution by a collective response of the bulk of the pro- 
duced particles to the shift in the center of mass in the initial state is not im- 
mediately consistent with a parton model approach. To achieve similar results, 
one must introduce "shadowing" and "anti-shadowing" of the bound-nucleon 
structure functions to achieve similar effects in the rapidity distributions. How 
this should preserve "wounded nucleon" scaling in a natural way is not clear 
from the perturbative approach, perhaps leaving room for simpler, heuristic 
models to point to issues where additional theoretical guidance will be needed 
in the future. 



As a final note: as this manuscript was being completed, the author became 
aware of similar approaches to understanding the RHIC d+Au data, based 
on "shifting" the bulk distributions in rapidity[15] or pseudorapidity[16,17]. 
It should be noted that there are basic differences. For example, Jeon et al 
[16] compare d+Au data with Au+Au data, applying scaling and shifts in 
dN/drj and shifts in pseudorapidity space, where the work here considers shifts 
in rapidity and no rescaling of the multiplicity beyond participant scaling 
between d+Au and p+p data. Still, it is intriguing to see similar ideas to those 
presented here appearing independently in both theoretical and experimental 
works. 

This work was supported under U.S. DOE grant DE-AC02-98CH10886. The 
author acknowledges insightful discussions with Mark Baker, Wit Busza and 
Rachid Nouicer. 
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